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ABSTRACT 


A generalized  measure  of  performance  is  defined  as  a weighted 
combination  of  the  ergodic  queue  length  distribution  where  the  weights 
are  general  functions  of  the  system  parameters.  The  paper  presents  a 
sequence  of  upper  and  lower  bounds  for  this  measure  of  performance  in  the 
G/M/l/N  queue  with  FIFO  discipline.  The  bounds  are  used  to  derive  a 
sequence  of  approximations  with  bounded  errors.  The  upper  and  lower 
bounds  are  shown  to  converge  to  their  corresponding  exact  values.  The 
technique  used  is  based  on  the  imbedded  Markov  chain  analysis  and  consid- 
ers only  subsets  of  the  steady  state  equations  to  derive  the  bounds. 
Initial  computational  experience  is  encouraging  and  has  indicated  that  the 
approximations  are  viable  for  heavy  and  medium  traffic  conditions. 


1 . INTRODUCTION 


Consider  a single  server  queue  with  a finite  capacity  N and  'first 
come,  first  served'  discipline.  The  arrival  epochs  form  a recurrent  proc- 
ess and  an  arriving  customer  leaves  the  system  without  service  if  the 
number  of  customers  in  the  system  is  N.  The  service  times  are  independ- 
ent and  identically  distributed  random  variables  with  negative  exponential 
distribution  with  mean  y . In  an  extended  Kendall  notation  a system  of 
this  type  can  be  identified  as  G/M/l/N.  This  queue  and  its  generalization 
G/M/s/N  (with  s servers)  are  combined  loss  and  delay  systems  occuring  in 
telephone  trunking  and  communication  problems. 

1.1  Preliminary  results: 

Let  t^,  t^,  ty  ...  be  the  arrival  epochs  in  the  system,  and 

U = t - t , be  independent  and  identically  distributed  random  variables 
n n n-1 

with 


P [U  < x]  = A(x) , x > 0 
n — — 

and  E[U  ] » X \ 
n 


Let  Q(t)  and  Qn  be  the  number  of  customers  in  the  system  at  time  t and 
just  prior  to  the  nth  arrival  respectively.  Because  of  the  exponential 
nature  of  service  times,  {Q^,  n = 1,2,3,...}  is  a finite  Markov  chain  imbed- 
ded in  the  stochastic  process  (Q(t),  t 0}  with  transition  probabilities 
P^ , 0<^i,  j <_  N,  where 


p[Vi " Uq„  - u 


J (i  + 1 - j)! 
0 


dA(x) , i + 1 > j > 1,  N>i>0, 


1 


clearly, 


iO 


i+1 
1 - I P 
2=1 


ij 


N > i > 0 


and 


P = P 
Nj  N-l , j 


N > j > 0 


Let  3/  = (q^,  q^,  — , q^)  denote  the  steady  state  probabilities  of  the 
number  of  customers  in  the  system  just  before  an  arrival  occurs.  Then  it 
is  well  known  that  ^ is  the  unique  solution  to  the  system  of  equations 


3/  I.  = 3/  ,3/1=1 


(1) 


where  P^  is  the  one  step  transition  probability  matrix  and  is  given  by: 


P = 


1 - a. 


1 - ao  - ai 


a„ 


N-l 

1 - 2 a.  a..  , 

n J N'1 


N-l 

1 “ q aj  aN-l 


al  a0 


al  a0 


where  we  have  written 


a 


i J 


^ e-yx  iM)/  dA(x) 


(2) 


0 


2 


* The  method  of  imbedded  Markov  chain  as  applied  to  G/M/l  is  due  to 

Kendall  [11],  and  further  use  of  the  approach  is  presented  in  Prabhu  [14], 
Keilson  [10],  Bhat  [2],  Laslett  [13]  and  most  recently  Raju  and  Bhat  [15]. 

1.2  A Generalized  Measure  of  Performance 

Measures  of  performance  are  functions  of  system  parameters  X,  y,  N, 
etc.  They  are  needed  in  operational  analysis  where  the  best  parameter 
values  which  optimize  system  performance  are  sought.  Measures  have  been 
used  since  Erlang's  work  "The  Rational  Determination  of  the  Number  of 
Circuits"  [3]  where  the  probability  of  loss  to  the  system  is  minimized. 

For  the  history  and  development  of  operational  models  and  their  measures 
see  Bhat  [1]  and  for  a recent  survey  see  Crabill  et  al  [5].  Typical 
measures  of  performance  are  expressed  in  terms  of  cost  of  service  per 
customer,  waiting  cost,  lost  customer  cost,  server  utilization,  etc. 
[Hillier  [8]].  Since  the  most  commonly  used  measures  can  be  written  as  a 
linear  combination  of  the  steady  state  probabilities,  a natural  extension 
will  be  the  generalized  measure: 

N _ 

ZN  = ^  I *  3i  qi  ’ ^ 

i=0 

I 

where  {a^}  are  general  functions  of  the  system  parameters.  For  applica- 
tions of  this  measure  see  Evans  [6]  and  Kotiah  [12] . 

This  investigation  considers  the  evaluation,  bounding  and  approxima- 
tion to  the  generalized  measure  as  defined  in  (3)  as  a necessary  step 
toward  optimizing  the  system  performance.  Early  experimentation  with  (3) 


3 


using  some  general  distributions  shows  that  ZN  is  nonconvex,  nonconcave, 

and  even  multimodal.  Since  optimizing  a generalized  function  (as  Z„)  is 

N 

likely  to  call  for  a large  number  of  functional  evaluations,  a fast  tech- 
nique to  evaluate  Z^  is  needed.  Moreover,  the  bounds  can  be  used  in  a 
Branch-and-Bound  framework,  and  the  approximations  can  be  utilized  to 
obtain  approximate  optimal  solutions.  Existing  analytical  techniques  are 
fairly  complicated  to  be  extended  for  operational  analysis,  but  a technique 
as  in  [15]  is  promising  for  the  fast  numerical  results  that  we  seek. 

The  paper  is  divided  into  five  sections.  Section  2 investigates 
the  ratio  properties  of  the  steady  state  distribution.  Such  properties 
will  facilitate  calculating  the  steady  state  distribution  and  reformulat- 
ing the  steady  state  equations  to  help  develop  the  bounds.  The  suggested 
sequence  of  upper  and  lower  bounds  to  Z^  is  presented  in  section  3.  A 
sequence  of  approximations  that  seems  to  work  best  for  heavy  and  medium 
traffic  is  given  in  section  4.  Computational  experience  for  the  approxi- 
mation is  given  in  section  5. 

It  should  be  clear  that  all  the  results  presented  here  using  the 
arrival  epoch  steady  state  distribution  can  be  easily  modified  for  the 
arbitrary  time  distribution.  This  is  true  by  virtue  of  Hokstad's  [9] 
result  which  shows  that 

qK  = y \-l  ’ K = 1’  2’  • ‘ • ’ N 

where  tt  is  the  arbitrary  time  steady  state  probability  that  there  are  K 
customers  in  the  system  [see  also  Takacs  [16]  and  Heyman  [7]]. 
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2.  RATIO  PROPERTIES  OF  THE  STEADY  STATE  SOLUTION 

In  this  section  we  focus  on  the  properties  of  the  ratio  of  steady 
state  probabilities 


ri  " qN-i+l/qN-i  ’ 1 1»  2,  ...  N. 


We  show  that  it  suffices  to  generate  a single  sequence  r^,  i = 1,  2,  ....  K, 
independent  of  N,  to  calculate  the  steady  state  distribution  for  any  of  the 
queues 

G/M/l/i,  G/M/l/i  +1,  ...  , G/M/l/N,  ...  with  the  same  general  inter- 
arrival and  exponential  service  time  distributions. 

N 

Theorem  1:  The  arrival  epoch  steady  state  distribution  {q^}^_p  of  the 
G/M/l/N  has  the  form 


q0  C0 

q±  = q^  , i = 1,  2,  ...,  N where  (4) 

ci  > 0 and  unique,  and 

N-i  N-i-1  K 

Ci  = “(Z*1  ‘ al  " “N-i  .n  Ci+j  " Z “ ,n  Ci+j}’  (5) 

j=l  J K=1  K+l  j=l  J 

CN  = “o^1  ~ ’ 

N N 

cn  = 1/{1  + E n c.} 

0 K-l  i=l  1 

Proof : The  proof  proceeds  by  induction  on  i,  i = N,  N-l,  ...»  1 where  we 
solve  the  steady  state  equations  (1)  recursively  and  then  apply  the 
normalizing  condition. 
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For  i = N,  we  solve  for  the  last  equation  in  £ = 5L.'  > 


% a0  + V-l  a0  " % 


Defining  cN  = a^/ (1  - a^)  we  get  qN  = cN  qN-1  and  the  theorem  is 
true  for  i = N.  Assuming  q^  = CK  % 1 for  K = N»  N - 1,  •••  » i + 1,  and 
solving  for  the  (i+l)st  equation  in  the  system  of  equations  <±'  = c^'  yields 


3.  ' 


°r  Vl  a0  + qi  “l  + + qN-l  Vi  + qN  Vi  = V N°ting  th3t 

K 

qi+K  = Ci+K  qi+K-l  = ^ 11  ci+i*qi’  and  ^arranging  we  get 

J-l 


N-i-1 


qi-l  “O  ~ qi^1  “ “l  ~ Z Vl  ,n  Ci+j  Vi  .n  Ci+j* 

K=1  3=1  J 3=1  J 


which  leads  to  the  form  given  in  (5)  and  therefore  the  induction  is 
complete. 


Now  applying  the  normalizing  condition  I q.  = 1 and  substituting 

i=0  X 


Ci  = Ci  c2  ...  ci  qQ,  we  get 


N K 

qn  = 1/ { 1 + t II  c,}  = c_. 
K-l  i=l 
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It.  is  known  that  {q.}  can  be  uniquely  determined  by  solving  equa- 
1 0 

tions  (1),  hence  their  respective  ratios  c = q /q  > K = 1,  2,  N 

is.  K K— 1 

are  also  unique. 


Define  r.  = cN_.+1  = qN_±+1/qN_±  , i > 2 


rl  CN  qN/qN-l  ’ 


r2  CN-1  ~ qN-l/,qN-2  ’ **• 


Then  from  theorem  1, 


i-1  i-2  K 

ri  = a0/[1  " ai  ai-l  ^ cN-i+j+l  " KfL  “k+I  cN-i+j+l^’ 


but  c„  r.  . , and  therefore 

N-i+j+1  i-j 


'i  • V11  - ai  - Vi  ”,  Ii-j  - Vi  'i-j1-  1 i 2 

j=l  K=1  j=l 


rl  " a0/(1  - V 


The  sequence  r^  can  be  generated  recursively  and  independently  of  the  system 
capacity  N.  That  is,  the  sequence  {r^}  for  G/M/l/N^  is  the  same  as  that  for 

G/M/l/N^,  G/M/l/N^ Hence  the  {r^}  is  a generating  sequence  for  the 

class  G/M/l/N  with  the  same  interarrival  and  service  time  distributions. 


To  study  the  behaviour  of  G/M/l/N^  for  different  system  capacities 
N , -C r i } provides  a computational  convenience.  It  suffices  to  generate  {r^} 
once  as  in  (7)  and  for  every  N _>  1,  {q.}N  is  given  by 
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distribution).  Note  that  from  this  table  steady  state  distributions  for 
systems  with  the  given  interarrival  and  service  time  distributions  and 


N <_  10  can  be  determined  based  on  the  following  lemma. 


COV 

.447* 

Table  1 

* 

.577 

1.000 

** 

2.134 

K=5 

K=3 

K=1 

V-1 

rl 

3.169 

3.684 

4 

5.049 

r2 

5.668 

5.252 

4 

3.118 

r3 

6.781 

5.984 

4 

2.094 

r4 

7.334 

6.328 

4 

1.784 

r5 

7.576 

6.487 

4 

1.707 

r6 

7.669 

6.556 

4 

1.689 

r7 

7.699 

6.585 

4 

1.685 

r8 

7.705 

6.597 

4 

• 

1.684 

r9 

7.706 

6.601 

4 

1.684 

ri  n 

7.706 

6.603 

4 

1.684 

^ 1 ^ tv  x _ 

For  Erlangian  family,  dA(x)  = - — r-r“ — j~ , cov  = 

(K  - 1)  • /K 

**  -P-.  Ax  -p  Ax 

For  Hyperexponential  family,  d A(x)  = 2p^  A e + 2p^  A e 


pl  ' 1 - p2  ' cov  - 1 


1 
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Lemma  1:  Let  {qf^}  and  {qf^}  be  the  steady  state  distributions  of  the 
two  systems  G/M/I/N^  and  G/M/l/N2  where  all  factors  other  than 
and  N2  are  the  same  for  the  two  systems.  Then  we  have 


q(1)  /qu)  = q(2)  /qN  ± , n - i,  n - i > o 

Nj-i+l  N^-i  N2~i+1  2 


(9) 


Proof:  Both  sides  of  (9)  are  equal  to  r^,  where  r^  is  independent  of  N^ 
and  N^. 


The  existence  of  the  generating  sequence  {r^}  can  also  be  observed 
from  the  results  of  Raju  [15] , where  a monotone  nondecreasing  sequence 
{a^}  is  used  to  calculate  the  steady  state  distribution  for  the  G/M/l/N 
system.  The  derivation  of  these  results  is  based  on  a fundamental  recur- 
sion used  in  inverting  the  system  of  equations  (1). 

Theorem  2:  [Raju  and  Bhat  [15]] 

The  steady  state  arrival  epoch  distribution  of  the  G/M/l/N  system  is 
given  by: 


_ aN-l-i  ~ aN-2-i 


1 - a0  a0 


q'i  ’ qN-l  aN_1  ’ qN 


^-1 


^-1 


where 


a0  = 1 


ax  = (1  - c^)/^ 


aK+l  * 3"  (aK  - Vl  “K-c1  • K - 1 

0 t=0 


According  to  Theorem  2 and  by  virtue  of  the  uniqueness  of  {r^},  r will 
take  the  alternative  form 


rl  1 - a. 


r a0(1  - V 

2 1 - a0  - Oj 


ai-2  ai-3  . . 

1 li3 


Often,  one  would  like  to  increment  N to  N + 1 or  decrement  N to  N - 1 to 
study  corresponding  effect  on  both  the  steady  state  distribution  and 
other  measures  of  performance.  Theorem  2 gives  rise  to  some  interesting 
properties  as  in  the  following  lemma. 

Lemma  2:  The  steady  state  distribution  for  G/M/l/N+1  system  is  given  by 
(N+l)  _ (N)  Vl 

"i  -Vi  17  • vi. 


(N+l)  . 


Proof : Using  Theorem  2, 


q(N+l)  = Vi  aN-l-i  _ aN-l-(i-l)  aN-2-(i-l)  aN-l 
1 aN  aN-l  aN 


(N)  Vl 


= a(N+1)  + T a(N+1)  - a(N+1)  + 

aN-l 

N+l 

y 

q0  + 1f1  qi  q0  + 

L» 

i-1 

- a(N+1)  + 

aN-l 

N 

y 

qo  + 

L i 

i=0 

,00  _ N+l  Vl 


hence  the  result. 
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3.  BOUNDS  FOR  THE  MEASURE  Z„ 

N 

Consider  the  optimization  model 


V m‘"  z»  ' Z ai  ’i  ’ a'  1 

q>0  1=0 


9'  E = 9 

9'  i = 1 


(12) 


For  a given  general  interarrival  time  distribution,  y,  and  N the  system 

of  equations  (12)  has  a linear  objective  function  and  linear  constraints 

N+l 

in  g.  The  feasible  region  is  a single  point  q*  e E which  is  the  solu- 
tion to  the  steady  state  equations.  It  is  known  from  the  theory  of  linear 
programming  [4]  that  a relaxation  of  the  feasible  region  yields  a solution 
which  is  a lower  bound  to  the  objective  function  Z^  of  (12).  A speciali- 
zation of  this  approach  is  used  in  Kotiah  (12],  to  prove  that  the  proba- 
bility of  zero  delay  for  M/D/2  is  larger  than  that  for  M/M/2  with  the 
same  traffic  intensity  and  to  obtain  bounds  on  the  average  queue  length 
for  an  M/M/ s mixed  queue.  Here  we  develop  bounds  on  the  generalized  measure 
using  a particular  relaxation  for  P . 

4.1- 

Definition:  Let  the  i relaxation  PR.  of  P„  be  defined  as 
i N 


N 

I q,  - 1 

i=0 


» q±  Z.  o 


Thus  the  i^  relaxation  reduces  the  constraint  set  to  the  last  i constraints 
along  with  the  normalizing  condition. 


Using  the  sequence  {r^}  and  Theorem  1 we  can  write  PR^  as 


PR . : min  = a q 
l N a 

q>0 

s.  t. 


•ri  qN-l  - % = ° 


r2  qN-2  " qN-l  U 


(14) 


ri  qN-i  ~ qN-i+l 


= 0 , 
N 

I q.  = 1 
i=0  1 


The  constraints  in  (14)  are  equality  constraints  and  can  be  substituted 
in  the  objective  function  to  give  the  concise  form  of  PR^. 


_ N-i-1  _ _ 

PR . : min  Zx  = I a . q . + q__  . h . 
i N . _ l N-i  l 

<j>0  i=0 

N-i-1 

I q . + q . h . = 1 
s.t.  . _ l MN-i  l 

1=0 


where 


(15) 


h.  = aM  . + aXT  . r.  + ...  + aM  r.r.  . ..r., 
i N-i  N-i+1  i N i l-l  1 


h.  = 1 + r.  + r.r.  . + ...  + r.r.  , •••r1, 
l ii  l-l  i i-1  1 


where  {r^}  are  known  for  the  class  G/M/l/N,  N 1.  Now  any  solution  of 


(15)  will  provide  a lower  bound  on  Z 


N 
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Definition:  Define  the  i^  lower  bound  LB.  of  Z_.  as  the  optimal  solution 

1 N 

to  (15).  Theorem  3 gives  the  explicit  form  of  the  lower  bounds. 


Theorem  3: 

>_  LB^  - min  (a^,  a^,  ...,  aN_i.  ]_*  N > i >_  1 

Proof:  PR^  is  a linear  programming  primal  problem  with  one  constraint. 

Defining  DR^  as  the  dual  problem  of  PR^  with  w as  the  dual  variable, 
we  have: 


DR^:  max  1.  w 
w 


s.t.  W < aQ 


w < 3l 


w < a . 

— N-i-1 

h . w < h . 

l—i 

By  simple  inspection,  the  optimal  feasible  solution  to  DR.  is 

w*  = min  {I0,  Ir  ...,  \_±_v  Vhi}’  hi  > 0 

* — * 

Also  since  w is  optimal  feasible  to  DR^,  an  optimal  solution  Z ^ 

to  the  primal  PR^  exists  (see  [4])  and 
* — * 

W Zi  — ZN 

Definition:  Define  the  i*"^1  upper  bound  UB^  of  Z^  as  the  optimal  solution 
to  the  following: 
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_ N-i-1  _ _ 

max  Z = I a . q . + q h . 

3>0  N i=0  1 1 N-1  1 


N-i-1 

s.t.  Z q.+q„.h.=l 
i-0  1 N-1  1 


Theorem  4 gives  the  form  of  the  upper  bounds , 


Theorem  4: 


ZN-  UBi  = Mx  {a0>  •••»  Vi-1’  hi/hi} 


Proof  of  this  theorem  is  similar  to  that  of  Theorem  3. 


Theorem  5 shows  that  the  sequence  LB^  i = 1,  2,  N 

ing  and  converges  to  ZN  for  i = N. 


Theorem  5: 


Z*  < Z*  < Z„  , 2 < i < N-l 

l-l  — l — N — — 


with 


ZN-ZN 


Proof : By  Theorem  3: 


zN  > Z* 


for  all  N-l  > i > 1 


Let  F.  be  the  set  of  all  feasible  solutions  to  PR., 
i i 

set  of  all  feasible  solutions  to  P„,  then  since  PR. 

N : 

to  PRi> 


F £ h £ Fi-i 


(16) 


is  monotone  incres- 


and  let  F be  the 
* is  a relaxation 
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with 


F = V 


hence 


and 


* * 

ZN  — ■ Zi  — Zi-1 


— * 

Z =*  Z 
N N 


A similar  result  can  be  obtained  for  the  convergence  of  the  sequence 

of  upper  bounds  UB.  to  Z„T. 

1 N 

Thus  Z^  is  shown  to  be  trapped  in  a sequence  of  intervals  that  can 
be  made  arbitrarily  small  by  virtue  of  the  convergence  property. 


min  A , < Z < max  A.,  i=  1,  2,  N 


0<J<  N-i 


j - N - 


0<j  <N-i  3 


(17) 


r a 0 1 J 1 N-i-1 

A.  = { _ 

J [ h^fu  j = N-i 

To  illustrate  the  applicability  of  the  bounds  consider  the  following 
profit  functions: 


(1)  The  coefficients  {aj}  are  monotone  increasing  in  j,  the 
number  of  customers  in  the  system,  e.g. , 


N 

Z1  = 2 jX  q.  , x > 0 

j=0  3 


15 


(2) 


{a^}  is  monotone  decreasing  sequence  in  j,  the  number  of 
customers  in  the  system,  e.g.. 


N 

Z = I (1  - jy)  q.  , y > 0 
j=0  J 


(3) 


is  the  sum  of  the  above  two  sequences. 


Z3  = 


Zl  + Z2 


Table  2 shows  the  bounds  to  Z , Z^  Z 3 for  E5/M/l/5  system  with  x = .4, 
y = .2  for  p = Ay  ^ = .8,  2,  4,  10.  The  generating  sequence  {r_^}  is  given 
for  each  p.  The  results  show  that  the  bounds  are  tighter  under  heavy 
traffic  conditions  since  they  rely  more  on  qN>  q , ...  rather  than 
3q»  •••  • They  also  show  that,  in  general,  we  should  expect  either  the 

upper  or  the  lower  bounds  to  be  useful. 

The  form  (17)  suggests  different  approximations  to  ZN»  For  example 
ZN  - 1/2 (LB ^ + UB^)  where  the  error  will  be  known  not  to  exceed  1/2 (UB.  -LB.) 
1 i £ N-l.  However,  an  approximation  of  a particular  importance  and  inter- 
pretation is  the  Z^  - A^__^  = h /Ik  as  shown  in  the  following  section. 
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Table  2 


p 

i 

r . 
i 

(LB.., 

/—S 

•H 

r3 

2 

(LB±. 

t 

'2 

, UE  ) 

z, 

(LB., 

UB  ) 

.8 

1 

.487 

(0.000, 

1.792) 

(.134, 

1.000) 

(1.000, 

1.950) 

2 

.639 

(0.000, 

1.668) 

(.270, 

1.000) 

(1.000, 

1.938) 

3 

.675 

(0.000, 

1.517) 

(.412, 

1.000) 

(1.000, 

1.930) 

4 

.681 

(0.000, 

1.317) 

(.562, 

1.000) 

(1.000, 

1.880) 

5 

.682 

(0.840, 

0.840) 

(.721, 

0.721) 

(1.561, 

1.561) 

2 

1 

1.638 

(0.000, 

1.839) 

(.075, 

1.000) 

(1.000, 

1.950) 

2 

2.430 

(0.000, 

1.800) 

(.119, 

1.000) 

(1.000, 

1.920) 

3 

2.769 

(0.000, 

1.778) 

(.141, 

1.000) 

(1.000, 

1.919) 

4 

2.891 

(0.000, 

1.765) 

(.152, 

1.000) 

(1.000, 

1.918) 

5 

2.929 

(1.756, 

1.756) 

(.156, 

0.156) 

(1.913, 

1.913) 

4 

1 

3.619 

(0.000, 

1.865) 

(.043, 

1.000) 

(1.000, 

1.950) 

2 

5.668 

(0.000, 

1.853) 

(.056, 

1.000) 

(1.000, 

1.920) 

3 

6.781 

(0.000, 

1.851) 

(.059, 

1.000) 

(1.000, 

1.910) 

4 

7.333 

(0.000, 

1.850) 

(.060, 

1.000) 

(1.000, 

1.910) 

5 

7.576 

(1.850, 

1.850) 

(.060, 

0.060) 

(1.910, 

1.910) 

10 

1 

9.607 

(0.000, 

1.884) 

(.018, 

1.000) 

(1.000, 

1.950) 

2 

15.603 

(0.000, 

1.882) 

(.021, 

1.000) 

(1.000, 

1.920) 

3 

19.421 

(0.000, 

1.882) 

(.021, 

1.000) 

(1.000, 

1.904) 

4 

21.810 

(0.000, 

1.882) 

(.021, 

1.000) 

(1.000, 

1.904) 

4 

23.233 

(1.882, 

1.882) 

(.021, 

0.021) 

(1.904, 

1.904) 
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A.  APPROXIMATIONS  TO  THE  GENERALIZED  MEASURE  Z, 


Here  we  develop  a sequence  of  approximations  for  Z^  based  on  A^ 
t h A 

Definition;  The  i approximation  Z^  (i=l,2,  N)  to  Z^  is  given  by 

the  solution  to  the  system 


min  Z a . q . 
|>0  N-i  J 3 


"l  %-l  - qN  = 0 


r2  %-2  ~ qN-l  = ° 


ri  qN-i  qN-i+l 


qN-i  + qN-i+l  + 


+ %=  1 


where  q^  are  the  approximated  steady  state  probabilities,  N _>  j _>  N - i. 

A 

Lemma  3:  Z^  = h^/h^  and  the  approximated  steady  state  probabilities 

A 

{ q j } are  given  by 


'"j 

• n 


q^N)  • d N - i < j < N 


otherwise . 


where  d = and  (qjN^)  is  the  exact  steady  state  solution 


to  the  corresponding  G/M/l/N  system. 
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Proof:  First  observe  that  equations  (18)  and  (19)  are  the  steady  state 
equations  for  the  corresponding  G/M/l/i  system  expressed  in  terms  of 


( i)  ^ 

ratios  (rj.  Let  (q  be  the  solution  to  G/M/l/i  system,  then 


r.  ,<«  - ,<» 


Hi-2  Hi-1 


= 0 


r.  q (i)  - q(i) 


"i  qo 


= 0 


hence 


~ (i) 

qi  = qj~ (N-i) 


will  satisfy  (18)  and  (19). 
Also  from  lemma  2 : 


qj 


, (N) 


(N-l) 

aN-2 

qj-l 

^-1 

N-i, 

q(±) 

. "i-l 

qj-(N-i) 

Vl 

a(i) 

- 

qj-(N-i) 

‘ qj 

= „ (N-2)  ^-3 

a. 


= qj~2 


N-l 


or 


d , where  d = 


^-1 


*i-l 


Also  system  of  equation  (18)  implies 


qN-i  “ ri  qN-i+l  ^ [ri  ri-lJ  1 %- 


i+2 
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(N-K)  *N-K-1 

J~'i‘  aN-l 


q^-K 


[ri  Vl  *•  ril 


-1 


J 


> 


Thus  the  quality  of  the  ith  approximation  (to  be  measured  quantitatively  in 
the  following  section)  depends  on  the  factor  d and  the  traffic  intensity 
p = A/y.  The  closer  the  d to  1 the  better  is  the  approximation.  A heavy 
traffic  situation  justifies  setting  the  first  few  probabilities  to  zero.  It 
has  been  observed  in  experimentation  that  d is  also  related  to  the  traffic 
intensity  such  that  heavier  traffic  provides  values  of  d which  are 
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closer  to  1.  The  following  table  shows  d for  different  levels  of  approxi- 
mations i = 1,  2,  3,  . . , N - 1 for  Erlangian  systems  (0  < cov  < 1),  and 


Hyperexponential  systems  (1  < cov  < °°)  under  various  traffic  loads. 


d 

X 

u 

cov 

N 

i=l 

i=2 

! 

i=3 

i=4 

i=5 

i=6 

i=7 

i=8 

i=9 

i=10 

2 

1 

.447 

10 

1 241 

1 074 

i nns 

1.008 

1.002 

1 000 

4 

1 

.447 

10 

1.044 

1.006 

1.000 

10 

1 

.447 

10 

1.006 

1.000 

2 

1 

1.000 

10 

1.332 

1.142 

1.006 

1.031 

1.015 

1.007 

— 

1.003 

1.001 

1.000 

1.000 

4 

1 

1.000 

10 

1.066 

1.015 

1.003 

1.001 

1.000 

► 

10 

1 

1.000 

10 

i m n 

i nnn 

2 

1 

2.134 

10 

1.416 

1.255 

1.173 

1.121 

1.084 

1.057 

1.036 

1.021 

1.009 

1.000 

4 

1 

2.134 

10 

1.112 

1.056 

1.031 

1.017 

1.010  1.005 

1.002 

1.001 

1.000 

1.000 

10 

1 

2.134 

10 

1 D9T 

i nnn 

X • uwu 

i 1 

j 

Using  the  bounds  (17)  it  can  be  easily  shown  that  the  approximation 

also  bounds  the  measure  Z,T  such  that 

N 


and 


Z.  > Z.T  if  Z,  > max  (a.), 

N 1~0<j<N-i-l  J 


Z . < ZXT  if  Z.  < min  (aj. 

1_  N 1 — 0<j<N-i-1  J 


The  following  lemma  shows  that  if  {a^}  is  a monotone  increasing  sequence 

the  approximation  will  bound  Z,T  from  above.  An  important  special  case  is 

N 

n th 

whenever  a^  = j , n=l,  2,  3,  ...  which  gives  the  n raw  moment  of  the 


distribution  of  customers  in  the  system. 


Lemma  4:  Let  {a..}  be  a monotone  increasing  sequence,  then 


ZN±Zi 


i < N - 1 
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Proof : 


N-i-1 


zn  = .z  aj  qj  = z aiqi  + 2 a q 

J=0  J j=0  2 2 j=N-i  2 2 


~ Vi-l  (1  ' z . qj}  + 1 ai  q< 

j=N-i 


j=N-i 


J 3 


’ ‘x-1-1  + j-ii  <aj ' 

- N _ _ 

- aN-i-l  + Z (ai  " aN-i-l)qi’  by  lernma  3 

j=N-i  2 w 1 1 2 

= Vi-1  " Vi-1  . I . qj  + h 

j=N-i  J 


= Z. 


since  I q.  = 1. 
j=N-i  2 
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5.  COMPUTATIONAL  EXPERIENCE 


The  approximation  procedure  as  presented  in  the  previous  section  has 

been  coded  and  used  to  obtain  computational  experience  for  E /M/l/N 

K 

(K-Erlangian)  systems  and  Hp/M/l/N  (Hyperexponential)  systems.  These  two 
families  of  queues  are  commonly  used  in  the  literature  since  they  cover 
the  full  range  of  the  coefficient  of  variation  (cov)  and  they  provide  a 
generating  formula  for  the  elements  ou,j  =0,  1,  ....  N.  This  speciali- 
zation does  not  simplify  the  approximation  procedure  but  does  simplify 
the  coding  effort. 


For  different  combinations  of  N,  X,  y,  and  cov,  a set  of  100  problems 
were  solved  for  the  i1"*1  approximation  i = 1,  2,  . N - 1.  The  cost 
coefficients  {a^}  were  randomly  generated  as  uniform  (0,  1)  variables  for 

A 

each  problem.  Thus  the  resulting  generalized  measures  ZN>  Z ^ are  on  a 
(0,  1)  scale.  The.  choice  of  the  different  parameters  was  designed  to  re- 
flect the  effect  of  traffic  intensity,  coefficient  of  variation,  and 
system  capacity  (N)  on  the  quality  of  the  approximation. 

Two  measures  were  used  to  evaluate  the  i*"^  approximation,  the  absolute 

A A 

deviation  | ZN  - Z^j , and  the  percent  deviation  | Z^  - Z | /Z  , i = 1,  2,  .., 

N - 1.  For  each  problem  set  and  each  i,  the  mean  and  variance  of  the 
above  two  measures  are  given  as 


DEV  (I) 
SD  (I) 
PDEV  (I) 
SP  (I) 


mean  abs.  deviation 
variance  abs.  deviation 
mean  percent  deviation 
variance  of  percent  deviation 
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Scaling  the  cost  coefficients  should  not  affect  the  percent  deviation  while 
it  should  scale  the  absolute  deviation  with  the  same  factor.  Tables  3,  4, 
and  5 present  our  computational  experience. 

Clearly  the  approximation  works  best  under  heavy  traffic  loads.  For 
p = 5,  the  2nd  approximation  (i  = 2)  gives  rise  to  a max  deviation  of  .0082 
and  a max  percent  deviation  of  2.5%;  for  p = 10,  max  deviation  drops  to 
.0019  and  max  percent  deviation  drops  to  .6%  over  the  range  of  cov  and 
N considered.  This  fact  combined  with  the  simple  formula  for  this 
approximation 


V-2  + r2  aN-l  + r2ri  3N 
1+r2+r2rl 


makes  it  an  efficient  approximation  for  heavy  traffic.  For  medium  traffic, 
similar  results  can  be  obtained  for  i = 3,  i = 4.  Over  the  range  of  cov, 
it  seems  that  the  approximation  performs  better  for  smaller  values  of 
coefficient  of  variation.  A comparison  of  table  6 with  the  appropriate 
tables  3,  4,  and  5,  shows  that  the  approximation  yields  better  results  for 
smaller  values  of  N. 


As  expected  from  the  derivation  the  approximation  can  not  handle  low 

traffic  conditions.  Clearly  an  approximation  that  uses  q^,  q^,  q^,  ... 

should  yield  good  results  in  which  case  the  i*"*1  approximation  will  be  a 

function  of  r , r„  , , ...,  r„  . rather  than  r.,  r„,  ...,  r.. 

N N-l  N-i  12  1 

Finally  we  believe  that  the  simplicity  of  the  formula  and  the  deriva- 
tion of  the  bounds,  the  approximations,  and  the  generating  ratios  for  the 
G/M/l/N  system  have  great  potential  in  the  optimization  of  this  class  of 
queues . 
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Table  3:  Hp/M/1/10  System 


* 

cov 


2.134 


2.134 


2.134 


2.134 


:cov  = 


1 


\ 

y 

i 

DEV ( I ) 

SDOIP 

PDEV(I) 

SP(I) 

1. 

l. 

2 

.1006 

.0055 

.2043 

.0245 

3 

.0790 

.0034 

.1582 

.0131 

4 

.0638 

.0020 

.1299 

.0089 

5 

.0516 

.0014 

.1055 

.0061 

6 

.0432 

.0011 

.0885 

.0050 

7 

.0305 

.0005 

.0618 

.0022 

8 

.0240 

.0004 

.0493 

.0016 

9 

.0175 

.0001 

.0352 

.0006 

10 

.0000 

.0000 

.0000 

.0000 

3. 

l. 

2 

.0209 

.0002 

.0504 

.0024 

3 

.0117 

.0001 

.0274 

.0008 

4 

.0073 

.0000 

.0180 

.0004 

5 

.0045 

.0000 

.0106 

.0002 

6 

.0029 

.0000 

.0070 

.0001 

7 

.0018 

.0000 

.0044 

.0000 

8 

.0010 

.0000 

.0025 

.0000 

9 

.0005 

.0000 

.0013 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

5. 

l. 

2 

.0082 

.0000 

.0246 

.0010 

3 

.0040 

.0000 

.0123 

.0003 

4 

.0022 

.0000 

.0067 

.0001 

5 

.0010 

.0000 

.0030 

.0000 

6 

.0005 

.0000 

.0016 

.0000 

7 

.0002 

.0000 

.0007 

.0000 

8 

.0001 

.0000 

.0002 

.0000 

9 

.0001 

.0000 

.0002 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

10. 

l. 

2 

.0017 

.0000 

.0057 

.0001 

3 

.0006 

.0000 

.0022 

.0000 

4 

.0002 

.0000 

.0008 

.0000 

5 

.0001 

.0000 

.0003 

.0000 

6 

.0000 

.0000 

.0001 

.0000 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

coefficient  of  variation. 
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Table  4: 

M/M/1/10 

System 

cov 

X 

y 

I 

DEV (7) 

SD  (I) 

PDEV(I) 

SP(I) 

1. 

l. 

l. 

2 

.1190 

.2410 

3 

.0898 

.0039 

.1820 

.0180 

4 

.0783 

.0033 

.1595 

.0154 

5 

.0651 

.0022 

.1333 

.0109 

6 

.0489 

.0014 

.1016 

.0071 

7 

.0443 

.0009 

.0917 

.0042 

8 

.0355 

.0006 

.0730 

.0028 

9 

.0231 

.0002 

.0476 

.0011 

10 

.0000 

.0000 

.0000 

.0000 

1. 

3. 

l. 

2 

.0088 

.0000 

.0230 

.0008 

3 

.0028 

.0000 

.0072 

.0002 

4 

.0010 

.0000 

.0027 

.0002 

5 

.0003 

.0000 

.0009 

.0001 

6 

.0001 

.0000 

.0003 

.0001 

7 

.0000 

.0000 

.0001 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

1. 

5. 

l. 

2 

.0020 

.0000 

.0061 

.0001 

3 

.0005 

.0000 

.0015 

.0000 

4 

.0001 

.0000 

.0003 

.0000 

5 

.0000 

.0000 

.0001 

.0000 

6 

.0000 

.0000 

.0000 

.0000 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

- 

10 

.0000 

.0000 

.0000 

.0000 

1. 

10. 

l. 

2 

.0003 

.0000 

.0011 

.0000 

3 

.0000 

.0000 

.0001 

.0000 

4 

.0000 

.0000 

.0000 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

6 

.0000 

.0000 

.0000 

.0000 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

Table  5:  E2/M/I/IO  System 


cov 

A 

y 

1 

DEV (I) 

SD(I) 

PDEV(I) 

SP(I) 

.707 

1. 

1. 

2 

.1232 

.0071 

.2512 

.0318 

3 

.0966 

.0054 

.1981 

.0239 

4 

.0776 

.0036 

.1583 

.0170 

5 

.0602 

.0023 

.1229 

.0102 

6 

.0510 

.0015 

.1041 

.0070 

7 

.0439 

.0009 

.0894 

.0045 

8 

.0342 

.0007 

.0708 

.0039 

9 

.0228 

.0002 

.0468 

.0009 

10 

.0000 

.0000 

.0000 

.0000 

.707 

3. 

1. 

2 

.0054 

.0000 

.0125 

.0004 

3 

.0015 

.0000 

.0034 

.0003 

4 

.0004 

.0000 

.0009 

.0002 

5 

.0001 

.0000 

.0002 

.0001 

6 

.0000 

.0000 

.0001 

.0001 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

.707 

5. 

1. 

2 

.0014 

.0000 

.0039 

.0000 

3 

.0002 

.0000 

.0005 

.0000 

4 

.0000 

.0000 

.0001 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

6 

.0000 

.0000 

.0000 

.0000 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

.707 

10. 

1. 

2 

.0002 

.0000 

.0007 

.0000 

3 

.0000 

.0000 

.0000 

.0000 

4 

.0000 

.0000 

.0000 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

6 

.0000 

.0000 

.0000 

.0000 

7 

.0000 

.0000 

.0000 

.0000 

8 

.0000 

.0000 

.0000 

.0000 

9 

.0000 

.0000 

.0000 

.0000 

10 

.0000 

.0000 

.0000 

.0000 

I 
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.707  3.  1.  2 .0060  .0000  .( 

3 .0013  .0000  .( 

4 .0003  .0000  .( 

5 .0000  .0000  .( 


.707 

5. 

1. 

2 

.0013 

.0000 

.0036 

.0000 

3 

.0002 

.0000 

.0006 

.0000 

4 

.0000 

.0000 

.0001 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

.707 

10. 

1. 

2 

.0002 

.0000 

.0005 

.0000 

3 

.0000 

.0000 

.0000 

.0000 

4 

.0000 

.0000 

.0000 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

1.000 

3. 

1. 

2 

.0083 

.0000 

.0208 

.0009 

3 

.0025 

.0000 

.0066 

.0002 

4 

.0009 

.0000 

.0023 

.0001 

5 

.0000 

.0000 

.0000 

.0000 

1.000 

5. 

1. 

2 

.0022 

.0000 

.0062 

.0001 

3 

.0004 

.0000 

.0012 

.0000 

4 

.0001 

.0000 

.0002 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

1.000 

10. 

1. 

2 

.0003 

.0000 

.0009 

.0000 

3 

.0000 

.0000 

.0001 

.0000 

4 

.0000 

.0000 

.0000 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

2.134 

3.  1.  2 

.0186 

M 9 

.0507 

3 

Bn » * B 

.0269 

4 

jjfl 

B ! t n 

.0141 

BBT!  5 

5 

B 1 1 n 

.0000 

.0000 

2.134 

5.  1.  2 

.0081 

.0000 

.0238 

.0009 

3 

.0037 

.0000 

.0111 

.0002 

4 

.0013 

.0000 

.0041 

.0000 

5 

.0000 

.0000 

.0000 

.0000 

2.134 

10.  1.  2 

.0019 

.0000 

.0069 

.0001 

3 

.0006 

.0000 

.0022 

.0000 

4 

.0002 

.0000 

.0005 

.0000 

5 

.0000 

.0000 

.0000 

.0000 
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